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I. INTRODUCTION 



Accurate path control of ships and underwater vehicles along prescribed 
geographical paths is a fundamental problem which is becoming increasingly important, 
particularly as the missions of ocean vehicles become more sophisticated with strict 
requirements for performance. In order for a control law to be able to perform its 
mission in a realistic operational scenario it has to be robust enough so that it can 
maintain stability and accuracy of operations in the presence of modeling errors and 
environmental uncertainties. The robustness properties of the design are particularly 
important due to the unpredictable nature of the ocean environment and the changes in 
the hydrodynamic characteristics of the vehicle during turning, changes in the forward 
speed, or operations in proximity to other objects in the area. For these reasons, there 
exists a need for the analysis of the robustness characteristics of a particular control law 
design and the establishment of a rational operational envelope based on stability and 
performance criteria. Previous studies (Parsons & Cuong, 1977) showed that gain 
adaptation is highly desirable due to changes in the linearized vehicle hydrodynamics 
with different operating conditions, such as depth under keel. The resulting adaptation 
scheme (Parsons & Cuong, 1980) required significant vehicle motion which may be 
undesirable when operating in restricted water, or in object recognition and localization 
tasks. Integral control techniques (Parsons & Cuong, 1981) proved quite effective, but 
neglected the behavior of the vehicle which becomes very important at low speeds and 
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hover operations. Model based compensators exhibit robust behavior under conditions 
of parameter uncertainty which is as good as the classical linear quadratic regulators for 
linear output feedback systems (Healey, 1992). Alternatively, sliding mode controllers 
exhibit very robust characteristics given an estimate of the parameter uncertainty and/or 
disturbances (Papoulias & Healey, 1992), (Yoeger & Slotine, 1985). Sliding mode 
control, however, does not offer an infinitely robust design, and it suffers from a series 
of bifurcation phenomena and loss of stability unless proper care is exercised (Papoulias, 
1991). 

In this work we analyze the problem of loss of stability of a path keeping control 
law under conditions far from nominal. We assume that the autopilot has been designed 
based on a nominal model, whereas the actual system is different. For demonstration 
purposes we employ a linear full state feedback control law, but the methods are quite 
general and can be used for other designs as well. The main loss of stability cases 
analyzed here occur in the form of generic bifurcations to periodic solutions 
(Guckenheimer & Holmes, 1983). We use center manifold reduction techniques and 
averaging in order to capture the stability properties of the resulting limit cycles (Chow 
& Mallet-Paret, 1977). Particular emphasis is placed on the control gains as the primary 
bifurcation parameters, since they are related to gain margins in linear control theory 
(Friedland, 1986). We make extensive use of numerical integrations in order to confirm 
the theoretical results. All computations in this work are conducted for the NPS 
autonomous underwater vehicle (Bahrke, 1992), and all results are presented in standard 
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dimensionless quantities with respect to vehicle length, 7.3 ft, and nominal forward 
speed, 2 ft/ sec. 
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H. PROBLEM FORMULATION 



A. INTRODUCTION 

The equations of motion of an ocean vehicle in the horizontal plane are presented 
in this chapter. A linear feedback control law is designed based on the linearized 
equations in yaw, sway, and rate of change of heading angle and lateral deviation error 
to provide path keeping. Loss of stability is examined for small changes from nominal 
in feedback gains, and system properties. 

B. EQUATIONS OF MOTION 

The mathematical modeling of a steering system of a vehicle consists of the 
nonlinear sway and yaw equations of motion. Newton’s equations of motion in a moving 
coordinate frame fixed at the ships geometrical center are 1 

m(v+r+x G f)=Y/+Y.v+Y r r+Y v v+Y 6 & S +Y.6 b 

s b 

stem 

- /c^(C)(v + Cr)|v + Cr|JC 

bow 

Ir+mxJv+r)=N-+N.v+Nr+Nv+N, 6 +N* 6. 

Z G v 'TV T V o s S 0^ b 

bow ( 2 ) 

- f C D h(0(v+Cr)\v+Cr\dC 

stem 



1 Equations are nondimensionalized with constant forward speed u,ship length L, and the dimensionless time 
being t.u/L. 
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where only nonzero terms are kept in the model and the symbols are explained in the 
nomenclature. The cross flow integral drag terms in the equations of motion are very 
important in low speeds, but for higher speeds their effect is much smaller and the 
steering response predominantly linear. Also for maximum maneuverability, the vehicle 
bow rudder is deflected at the same amount and opposite to the stem rudder. 

Equations (1), and (2) can be written in the form 



v=a n v+a l2 r+b 1 b 


(3) 


f=a 2i v+a 22 r+b 2 6 


(4) 



where 

Da n = (I z -N.)Y v -(mx G -Y.)N v , 

Da n= (A -fyX* - Y r ) ~{mx G -y r ){n r -mx c ) , 

Da 21 = (m-YJN v -(mx G -N,)Y v , 

Da 22 = (m-YJiN'-mxJ -(mx c -NJ(Y r -m) t (5) 

M>i = » 

Db 2 = ( m -YJ(N,-N h ) , 

D = (I z -N,)(m-Y i )-(mx G -Y.)(N r -mx G ) . 

The transfer function between rudder angle 5 and yaw angular velocity r is obtained from 
Equations (3) and (4). 

r_ _ b 2 S+(a 21 b l -a n b 2 ) ^ 

6 S 2 -(a 11 +a 22 )S+(a n a 22 -a 12 a 21 ) 

For low frequency maneuvering motions this second order equation can be 
approximated by expanding in Taylor series and keeping only the first order terms only 
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or r = ar+bb 



(7) 



r_ _ b 
6 s-a 

where 



( Q ll Q 22~ Q 12 a 2l)( Q+21 ^ ~ a i 

( a il +a 22)( a 21^1 -a i "A A I a 22 ~ a i2 fl 21^ 
AA^lA)" 

(a n +a 22 )(a 21 b 1 -a n b l -<3 n & 2 ) +b 2^ a n a 22 ~ a n a 2i) 



( 8 ) 



Equation (7), Nomoto’s Equation is very useful in control system design since no 
sway velocity feedback necessary, and it gives the fundamental turning performance. At 
steady state, Equation (6) turns out to be 



f = 0 

ar+bb = 0 ^ 

b x 

r = — o 
a 

From Equation (9) it can be seen that the relationship between yaw rate and rudder 
control angle are linear with a slope -b/a. Experimental results show that the actual 
relationship is not linear. Large increases in rudder angle fail to increase the yaw rate 
according to the amount predicted from Equation (9). 

The linear change in the time derivative of yaw rate due to rudder angle can be 
augmented by a nonlinear term of the form , a 3 r 3 ’. This is introduced to model the 
appropriate speed loss during turning. The term ’a 3 ’ has the same sign as ’a’. The new 
yaw equation is therefore 
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( 10 ) 



r=ar+a 3 r 3 +b6 




Figure 1 Difference between Nomoto’s equations and experimental results. 



The model becomes complete with the expression of yaw rate and the inertial 



deviation rate from the commanded path. Therefore the coupled nonlinear equations of 
motion for the marine vehicle in the horizontal plane are 


= r 


(11) 


r = ar+a 3 r 3 +bb 


(12) 


y = sinijj 


(13) 



This system of equations forms the basis for the controller design. 
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C. CONTROL LAW 

Equations (11), (12), (13) govern the steering control of the model used in this 
section. The control law can be expressed as, 



6 



a^tanhA 

6 sat 



where in the vicinity of ^ = r = y = 0we have 



( 14 ) 
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6 0 = K^ + K r r + K y y (15) 

5 is the rudder angle and K* ,K r , K y are the control gains of the system. The linear 
control law is 6 0 . The rudder angle 5 is defined by a hyperbolic tangent function to 
include the saturation to our problem as shown in Figure (3). Saturation occurs at 5 sat 
which is the saturation limit generally can be taken as 0.4 rad. The linearized form of 
equation of motions in the vicinity of ^ = r = y= 0 are 



ijr = r 


(16) 


r =ar+bb 0 


(17) 



y = ^ 



(18) 




Figure 3 

These equations can be expressed in space state form as 

X = XX 
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A is the Jacobian matrix of the system and X is the state vector 



A = 



0 1 0 
bK^ a+bK r bK y 

1 0 0 



( 20 ) 



X = [ *, r, y ] r 



The characteristic equation of the matrix A is 



( 21 ) 



X 3 -( a +bK r )X 2 -bKy X -bK y = 0 



( 22 ) 



The desired characteristic equation has the form 



A. 3 -*-a 2 A 2 +a 1 X+a 0 =0 (23) 

The gains can be generated in terms of the coefficient of actual and desired characteristic 
equation from Equations (22) and (23), 




( 24 ) 



K 



a 2 +a 

b 



( 25 ) 




( 26 ) 



The desired characteristic equation can be written with respect to the desired 
natural frequency and some optimum coefficients. The ITAE criterion for a third order 
equation is (Dorf, 1992) 
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S 3 + 1.75 u n s 2 + 2.15 o*S + <*l =0 



(27) 



Therefore the control gains can be calculated for a given natural frequency, as 



a o = <■>» 

= 2.15 w* 
cc 2 = 1.75 o) n 



(28) 



D. LOSS OF STABILITY 

1. Introduction 

The control law guarantees stability if all real parts of eigenvalues of the 
jacobian matrix are negative. Equation (23) is stable with the chosen control gains. A 
small perturbation in the control gains or coefficients of the system affects its stability. 
In this section we presented the computations for the critical value for change in the 
gains and the coefficients of the system. 

2. Stability Considerations 

For a system with a third order characteristic equation, the Routh Hurwitz 
criterion (Dorf, 1992) requires for stability 

a 2 a^a 0 (29) 

a. Perturbation in 

The change in is defined by a coefficient c, where c is any real 
number. The linearized equations of the system for a small change in are 



11 



\jf =r 
r =ar+bb 0 

y = y\f 



where the control law is 



6 o = cK ^ +K r r+K y y 
Therefore the characteristic equation is 



\ 3 -(a + bK r )X 2 -cbKyX-bK y = 0 



Equation (32) should meet the Routh Hurwitz criterion for stability. 



(a + bK r )cbK^-bK y 

The critical value of c where a stability change occurs in the system 
Equation (34) 



Using the coefficients, the critical c is 

C- = 0.2658 

cntr 



(30) 

(31) 

(32) 

(33) 

is presented in 

(34) 

(35) 
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Figure 4 c crit vs natural frequency for K^. 
b. Perturbation In K, 

Similarly for perturbation in K r the characteristic equation is 



~(a+cbK r )\ 2 -bK^X -bK y = 0 



(36) 



The critical value of c, 



a 0 +a 1 ^ 



cm, 



a,(a 2 +a) 



(37) 



or 



co +2. 15a 

c . = 

enr *' 2.15 ( 1.75 u n +tf) 



(38) 
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Ccritical 



20 




wn 

Figure 5 c crit versus natural frequency for K^. 



c. Perturbation In Ky 

The results for changes in K y are presented similarly.The characteristic 



equation, is 



\ 3 ~(a + bK r )\ 2 -bKy\ -bcK y = 0 



( 39 ) 



The critical value of c is 



a i tt 2 



cnl. 



( 40 ) 



y a. 



c„, = 3.7625 



cnt 



( 41 ) 
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Ccritical 
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Figure 6 c crit versus natural frequency for K r 
d. Perturbation in a 

After similar calculations the critical value of c is 



(a 2 +a) ctj -a 0 



cnt„ 



at t. 



(42) 



1.28484 a) 



C crit„ 



^ + 1 



(43) 
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Figure 7 c crit versus natural frequency for a 
e. Perturbation in b 

After similar calculations as with the other cases the characteristic 



equation is, 

A 3 + (c(a 2 +a) -a)a 2 +ca 1 X +ca Q =0 
and the critical values are 



(44) 



or 



cm. 



cc 0 +a a 1 
a x (a 2 +a) 



(45) 
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( 46 ) 



oi n +2.15 a 

Ccril > " 2.15(1.75 G) n +a) 

Examining Equation (46), we note that all the coefficients of polynomial 
must have the same sign, if all the roots are in the left hand plane. The necessary 
conditions for stability are, 



c^O 

c(a 2 +a) -a>0 



(47) 



The system is unstable for all values of the natural frequency if c is less then zero. When 
the natural frequency is in the interval 0 < w n < | a/ 1.75 1 stability changes from stable 
to unstable with an increase in c (Figure (8)) for positive values of c. In the interval 
a) n > 1 2.15a | stability changes from unstable to stable while c crosses the critical point 
again for positive values of c. In both cases a pair of complex conjugate eigenvalues 
crosses the imaginary axis. Stability is also lost when c becomes negative, but this is 
associated with a real eigenvalue crossing zero. In this work we concentrate on the c > 0 
case since in applications it is unlikely that a change in the sign of b will occur. 
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Ccritical 




Figure 8 c crit versus natural frequency for b.The c crit goes to infinity asymptotically at 
« a =-a/1.75. 
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m. HOPF BIFURCATION 



A. INTRODUCTION 

Hopf bifurcation is the simplest bifurcation, in which under the variation of a single 
control parameter a stable focus equilibrium bifurcates into an unstable focus surrounded 
by a growing limit cycle. As c crosses the critical value one pair of complex conjugate 
eigenvalues of the linear system matrix crosses transversely the imaginary axis. In this 
generic Poincare’ - Andronov - Hopf bifurcation a family of periodic solutions coexists 
with the stable / unstable nominal equilibrium state. Locally as c approaches c crit the 
periodic solutions are located on the two dimensional Euclidean plane spanned by the 
eigenvectors of Jacobian matrix of the system which corresponds to the critical pair of 
eigenvalues. In this chapter stability 
properties of the periodic solutions are 
established. In order to establish those 
properties the main nonlinear terms that 
dominate the system are isolated. Center 
manifold theory is used to reduce the flow 
to a two dimensional manifold. The 
method of averaging is applied to the 
reduced system. Hopf bifurcations are examined for small changes in each gain and 
system dynamics coefficients. 



mupnary 
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B. THIRD ORDER EXPANSION OF THE SYSTEM EQUATIONS 

1. Perturbation in K* 

In this case the equations of motions are 



where 



♦ = r 
r = a 
y = sinifj 



r = ar+a 3 r*+bb 



( 48 ) 



6 = 6_tanh 

sat 



^ sat ) 



( 49 ) 



S 0 = cK^+K, r +K y y 

The trivial equilibrium state is characterized by i/'=r=y=0 .Taylor expansion of the 
nonlinear terms about the equilibrium gives 



sinijr 

6 



+0(5) 

3! 



So- 



38 



■0(5) 



sat 



( 50 ) 



The system equations are 



i|f = r 

r = ar +a 3 r 3 +bb Q - 



b 8; 



36 



sat 



I v 

y = ^-- 7 - 

O 



( 51 ) 



20 



In state space form they are written as 



X =AX+g\X) 



( 52 ) 



where A is the linearized system matrix, and g 3 (x) contains the nonlinear terms, 



A = 



0 



0 



cbK^ a + bK r bK y 



0 



0 



( 53 ) 



gfc) 

g!(x) 



g s(*) 



= 0 



b >. 3 
°0 



36 

l! 

6 



sal 



where the S 0 3 term is computed as 



( 54 ) 



6 0 = cK^ 3 +K r r* +K y y 

3K^y 2 (cK^ + Kr) + 6cK v K r K y ^ry 



3c 2 Kl(K r r + K y y) + 3K?r 2 (cK^+K y y) 



( 55 ) 



With the computed gains from equations (24), (25), (26), Jacobian matrix is 



A = 



0 1 0 
-ca, -a 2 -a 0 

1 0 0 



( 56 ) 



The eigenvalues of matrix A are computed at the bifurcation point c crit from equation 
(34) where a pair of complex conjugate roots with real parts are obtained. 
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( 57 ) 



*i 



\ 



<£r 



a„ 



K = -j 

= -a 



a,. 



_o 

N “2 



3 “ 2 

For the above system of equations a transformation matrix of eigenvectors can be 
introduced, 



The transformation 



T = 



0 



0 - 



\ 



a. 



o, 



-a. 



0 



\ 



a, 

— -1/a, 

a o 



x-Tz , z=T~ l x 

transforms the system into its normal coordinate form, 



z = T- l ATz + T' l g\Tz ) 



where, 



( 58 ) 



( 59 ) 



( 60 ) 
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T -i = 



(olA 3 ! 2 



\ a oJ 



/ 3 \ 
( a 2 + a 0 ) 



( a n \ 3 ' 2 



a. 



v v 



ou 



a 2 +a 



o 



\ «o \ 



ttr, 



a-, 



0 



a. 



a- 






a 



2 

a O > 



a r 






(61) 



and 



rUr = 



o 



N 



a- 



0 



a. 



o 

N a 2 



0 



0 



0 



0 



CC n 



(62) 



The transformation of physical variables to normal coordinates is 

x= T z , 



^ = z 1 +z 3 



r = - 



a o 

^2 a 2^3 

a 2 



(63) 



y = 






a 2 1 



a 



ct n 
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1 



0 



1 



Y 

r 

. y 



o 




o 




-a 



2 




*1 



*2 

Z 3 



( 64 ) 



The coordinate z 3 corresponds to the eigenvalue X 3 which is real and negative. 
X,, X 2 denote the complex conjugate pair with zero real parts. At the critical point two 
distinct eigenspaces are spanned by the eigenvectors associated with the two sets of 
eigenvalues. These eigenspaces can be viewed as local approximations of the invariant 
manifolds the center and the stable manifold. The system exhibits its essential 
bifiircational behavior on the centre manifold associated with X, and X 2 . Centre manifold 
theory reduces the flow to a two dimensional manifold (Hopf bifurcation). According to 
centre manifold theory the coordinate z 3 is expressed in terms of z, and Zj, but this 
expression is higher order. Equation (64) can be written as, 



t = Zi 



r = - 



\ 



a 



a„ 



*2 



y = 






\ “o 



( 65 ) 



Substituting the equations (65) to (60) the system equations in normal coordinates 
become 
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( 66 ) 



*t = " 



Ctr 



- z 2 +r u zf +r u ziz 2 + r uZl zt+r u Z2 
\ a 2 



Z2 = 



\ a 



a 0 32 23 

— Z,+ r 21 Z 1 +/- 22 Z 1 Z 2 +/- 23 Z 1 Z 2 + ^*24 Z 2 



The terms r y are computed later from our equations at the bifurcation point c crit . For 
values of c in the neighborhood of the critical point, equation (66) becomes, 





/ 


\ 


II 

TO 

^1 


0) ; G + 


«0 


\ 


“2, 



*2 = 





\ 


/ 


a o 


Cl) € + 






a 2) 



z 2 +r u z 1 3 +r 12 z, 2 z 2 + r n z x z\ +r, 4 z 2 



Zl + V ' eZ 2 + r 21 Z l +r 22 Z l Z 2 +r 22 Z l Z 2 +r 24 Z 2 



(67) 



The parameter e is the difference between c crit and c or, c=c crit +e. The terms /?' and 
are the derivative of the real and imaginary parts of eigenvalues with respect to c 
evaluated at c=c crit , and computed from a perturbation series approach. The perturbation 
is expressed as follows from equation (34), 



a o 

c= +e 

a i a 2 



( 68 ) 



The characteristic equation can be written as follows for a small change for in the 
neighborhood of c crit . 



a 

X 3 + a 2 X 2 + (— +cCje)X +a 0 = 0 



(69) 
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The solutions of equation (69) are presented in equation (70) 



^1,2 “ b i e+ J 



X 3 = -a 2 +a l e 



/ 


— \ 




a o 


A 


+^6 

a 2 j 



( 70 ) 



The variables and Cj can be found by substituting equation (70) to (69) and neglecting 
the successive powers e 2 ,e 3 . 



\ = - 



a 2 a l 



2( a 0 + a 2 > 



U 2 a U 



C l = 



(71) 



CL r 



0 

M a 2 



2(a 0 +ap 

(S' and co' can be expressed by the definition of the derivative where e is zero in the limit, 



/ b.e- 0 
p 7 = — = b 



GO = 





e 


( 






«0 




«2 



+ Cj€ 



d r 



(72) 



N «2 



Equation (72) can be written in terms of the desired namral frequency, 



p 7 = -0.6841 a) n 
co 7 = 0.5171 «_ 



(73) 
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a. Calculations of r u Terms 

The terms are generated by the last term of equation (60), where 
T‘g 3 (Tz) is the third order nonlinear part. From equations (61) and (54) 



T- l g\Tz)= 









a 2 /2 ( a 2 +<£ o) 


a 2^ 2 ( a 2 + a 0 ) 


0 


3/2 


a 2 


0 



h. 

a o 

j_ 

a 2 

a. 



a 0 *2 



\3/2 



3 b *3 
Z 2 



1 3 

t 1 



( 74 ) 



where, 



= ±i- 



a oV 3 , 
— Zi +3 

V a 2 



a r \ 5 ' 2 



ou 
V 2) 



2 0 
az x z 2-3 



/ a \ 2 

a o 

a. 

V u 2 



2 2 3 

a z x z 2 + a 



<X A \ 3/2 



a. 
2) 



zl) ( 75 ) 



After expanding the first row of equation (74) and organizing the terms, we get 



a 2 3 

r n = ( a 2 + a o) 



a- \ 2 ^ 






( 76 ) 



/a r V' 2 



r i2 = - 



a. 

V 2J 



, 3 v a 

(a 2 + a 0 )— — 

b 2 b 



( 77 ) 



sat 



, 3 x a 
r 13 = (« 2 + a o) 



2s 2 » 



( 78 ) 
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(79) 



a 



r !4 = - 



— («2 + s )( fl 3 + ; 2 ) 

a oJ 3* 2 «L 



Similarly after expanding the second row of equation (74) and organizing, 



r 2i = -(a 2 + S) { 



3/2 5/2 

3 \ t 1 «0 1 a 2 } 



3* 2 6Laf 






0/3 . a 

cu ft 5 



sat 



a 



r 23 



1/2 
0 , 3 



a 



3/2 



(a 2 + a 0 ) 



a 






r 24 " 



«2 + a 0/ a- 



a. 



(a,+ 2^~ ) 
2 X" 5 

sat 



3b 2 6 



we get 

(80) 

(81) 

(82) 

(83) 
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b. Averaging 

Equation (67) can be written in a simpler form , 



Zj = P'ezj- 



*2 = 





\ 


/ 


«0 


G) e + 




N 





/ 


— \ 


/ 


0 


Cl) € + 




N 


“2) 



Z 2 +F 1 (z 1 ,Z 2 ) 
Zj+ P'eZj +F 2 (z,,Z2) 



( 84 ) 



where F, and F 2 contain the third order expansion terms. The use of polar coordinates 
makes it possible to decouple equation (84) . We use the polar transformation 
z,= R cos(0) , z,=R sin(0) 

After some algebra, equation (83) can be expressed in terms of R and 6, 

R = p / ef?+F 1 (f?cos(6),i?sm(0))cos(0) +F 2 (/?cos(0),f?sin(0))sm(0) (85) 



F0 = 





( 


— \ 






«0 / 






— + ( j ) e 




In 


a 2 J 



(86) 



R +F 2 (i?cos(0),/?sin(0))cos(0) -Fj(i?cos(0),Fsin(0))sin(0) 

Equation (85) can be written as follows, 

R = p 'eR + P(B)R 3 

where P(0) is 2ir periodic in the angular coordinate 6, 



(87) 



P(0) = r n cos 4 (0) +r 12 cos 3 (0)sin(0) +r 13 cos 2 (0)sin 2 (0) +r 14 cos(0)sin 3 (0) 
+r 21 cos 3 (0)sin(0) +r 2 ,cos 2 (0)sin 2 (0) +r 23 cos(0)sin 3 (0) +r 24 sin 4 (0) 



( 88 ) 



Equation (87) is averaged over one cycle to obtain an equation with constant coefficients. 
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R = $ f eR+KR 2 



(89) 



where K is defined by, 



2n 

K = -LJ* F(0) dd 



(90) 



Equation (90) is simplified as equation (91) after evaluation of the integral. 



^ g [^ r n +r i3 +r 22 + ^ r 24 1 



(91) 



Equation (91) can be expressed in terms of the natural frequency, 



K = 0.452 [3a 3 o>„ + 



u » , 3 a 1-75 2 , 

(to + q +<z ) col 

1.75 * 2 " 



50 / 



(92) 



Equation (87) can be expressed as, 



0 = 






a 



— + o )'e+F(G)R 2 



a~ 



(93) 



where, 



F(0) = r 2l cos 4 (0) + r 22 cos 3 (0) sin(0) + r 23 cos 2 (0) sin 2 (0) + r u cos(0) sin 3 (0) 
-r u cos 3 (0)sin(0) -r 12 cos 2 (0)sin 2 (0) -r 13 cos(0)sin 3 (0) -r 14 sin 4 (0) 



After averaging equation (93), a constant coefficient M is formed, 
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(95) 



e = 



OC r 



a. 



+ ( j) / e +MR 2 



where M is defined by, 




(96) 



or in terms of natural frequency, 



M = 0.7949 co* [ 3 ~ — - 2 0256 

( 1.75a co„ 



1 



> z b 2 



sat 



to. 



2 



)( 



CO, 



1.75 1.75 



-a 2 )] 



(97) 



2. Perturbation in K r 

The calculation of the formulas for the parameters K r , K y , a and b is similar 
to the calculation for K^. Therefore in the following sections only the main results are 
presented. The system equations for a change in K r are 

iJj = r 

f = ar+a 3 r 3 +bb (98) 

y = simir 



The Jacobian matrix is 



6 = 



6 ,tanh 

sat 



b„\ 



\^satj 



6 o = W+cKs+Kj 



(99) 
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A = 



0 1 0 

-a, [(1 -c)a-ca 2 ] -a 0 

1 0 0 



( 100 ) 



The eigenvalues of the Jacobian matrix are calculated at the bifurcation point from 
equation ( 38 ), 



* 1,2 




( 101 ) 



The state vector is presented in normal coordinates where Z3 corresponds to X 3 , 



V = z x 

r = ->[«X *2 




The system equation in normal coordinates are, 



( 102 ) 



*1 = % +r u z , 3 +r„z 1 2 z 2 + + 

h = (« 7 e+^) z 1+ (J'e^ + r 21 z? + r 22 z i Z2 + r 23 Z 1 z\ + r^z 2 



( 103 ) 



where, 
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(104) 



P 7 = - 



to' = 



a 1 (a 2 +a) 

2(a]+al) 
a 0 a\ /2 (a 2 +a) 

2(a\+ctl) 



= -0.45428 (1.75 c o n +a) 



= -0.1441(1.75co +a) 



The coefficients of the third order expansion terms are, 



r n = 



a o a i 



a 0 + a ! 



3b 2 5 



-a, - — 



sax 



r n = 



a 



a o a i 



7/2 



l2s2 2 

b &sat a 0 + a 



r i3 



a 



a o a i 



b 2 6 



2 2 3 

sat a 0 + a l 



5/2 



r !4 = 



«o a i 



a 0 + a l 



a 3 






3^ 2 C 



r 21 = 



a 0 + «i 



a 



11/2 



2 1/2 
a 0 a l 



3^ 2 6L 



r 22 



a 



a. 



£ 2 8L a n + a 



2 

' sat ^0 



(105) 



(106) 



(107) 



(108) 



(109) 



( 110 ) 
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(Ill) 



r 23 = 



a 



a 



9/2 



l2 S 2 2 3 

* 5 »r a 0 + aj 



a, 



r 24 



a 0 + a 1 



a 3 + 






The cubic coefficients K and M are obtained after averaging , 



ls:=0.052824[13.8675a 3 co^-0.5(o n +^H( (O 3 + 2.i5aa)J)( 2 .15 

2 j2 



b 2 d. 



sat 



M = 0.01 142 [20.3337 a 3 o)^ -0.7331 o> n 

20.9727 w* , 

+ -(a -4.6225 o) n )(2.15(o a +a 

* 2 sL 



3. Perturbation in Ky 

The formulas for perturbations in K y are presented in 
equations of motion are, 



i|r = r 

f = ar+a 3 r $ +bS 
y = sini|r 



5 



= 8 



sat 



tanh 



6 r 






\ 



/ 



6 0 = +K r r+cK y y 



The Jacobian matrix, 



( 112 ) 



u 2 n + a 2 )] ( 113 ) 



( 114 ) 

)] 



this section. The 



( 115 ) 



( 116 ) 
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0 

■a 

1 



1 0 

-«2 ~ ca 0 
0 0 



The eigenvalues at the bifurcation point , 



( 117 ) 



^ 1,2 " + j \[*\ 

X 2 = a 2 

The state equations in normal coordinates, 



'I' = Zi 

r = 'fa *2 




The normalized equations of motion, 

- = V'ez^v'e+Jal) % *r n z\ + r 12 z 1 2 z 2 + r 13 z,z 2 +r 14 z 2 

^ = (u'e+Ja^) z,+ p / ez 2 + ^ 21 Zi 3 + r 22 z 1 2 z 2 + r 23 z 1 z 2 

where /?' and co’ are given by, 



(118) 



(119) 



(120) 



P' = 



O) = 



dr 



2(cCj + cc 2 ) 



= 0.0959 to 
2 ■ n 



a 0 a 2 



2 /a^(a 1 + a 2 ) 



= 0.1145 o) 

2s " 



The coefficients for the third order terms, 



( 121 ) 
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( 122 ) 



r n - 



“l a 2 



2 2 
aj + a 2 



1 “i 



6 3 * 2 6 L/ 



r !2 = - 



5/2 

a «i “ 2 

^ 25 L “l +a 2 



(123) 



r !3 = 



fl 2 « 1« 2 
b 2 b] at ctj + a 2 



(124) 



r !4 = 



3/2 

a l a 2 



a i +a 2 



^3 + ' 






3 b 2 6 



2 X 2 



saf 



(125) 



r 21 = 



a 2V a i 



a 2 + otj 



1 “l + a 2 



3 ft 2 8 L a 2 6 



(126) 



r 22 = 



a 



CCi 



b a x + a 2 



(127) 



r 23 ” 



a 



a 



5/2 



b 2 v sa t “l + a 2 



(128) 



a, 



r 24 



aj + a 2 



a~ 



3 * 2 ««r 






(129) 



The cubic coefficients K and M , 
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2.15 cd 



K = 0.05155 [ 6.45 <z 3 cd„ -0.875 g>„ + (1.75 co n +a)(2.15 cd„ + a 2 ) 



2 b 2 



sat 



M = 0.3969 <z 3 cd* 



^2^-o>„(1.75ao>,-2.15 U „)(2.15G>^a 2 ) 

*> 2 «L 



4. Perturbation in a 

The equations of motion, 

r 

car+a 3 r 2 +bb 

sint|i 

where, 




The Jacobian matrix, 



6 = 6 ,tanh 

sat 



V 



b 0 = K^+K r r + K y y 



0 1 0 

A = ! - ce, fl(c-l)-a 2 -a 0 
1 0 0 



The eigenvalues at the bifurcation point, 




* J \j a i 

ct n 



CL 



and the state equations in normalized coordinates, 



(130) 



(131) 



(132) 



(133) 



(134) 



(135) 
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(136) 



t = z x 
r = -/«“ ^ 




The normalized equations of motion in the neighborhood of the bifurcation point, 



*1 = P /ez 1 - ( 0>/e+ /“l ) h. +r n z \ +r n z lh + r i3 Z l Z 2 +r u Z 2 

h = (o'e+^oQ z !+ p'ezj +r 21 z 1 3 +r 22 z 1 2 z 2 + r 23 z 1 z 2 2 +r M z% 



(137) 



p ; = -- 



G) 



a, 



a t + a 0 
3/2 

a o a i 



aj + a 0 



= -0.454a 



= -0.1441a 



(138) 



The coefficients of nonlinear terms, 



r n = 



«o a i 



3/2 



a 0 + 



5/2 

a, - 



3^ 2 C 



(139) 



a n a 



7/2 



r i2 = 



O^l 



« 2 s 2 2 3 

a 0 + a l 



(V^l( a 2 +a )~-~=) 
a, 



(140) 



— — ~r ( - (a 2 + a ) 2 a 2 - a 2 0 + 2 a , a 0 ( a 2 + a) ) (141) 

*> 2 $L«o + «i 
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1 



r i4 = 



3/2 

«o a i 
2 3 

a 0 +a, 



{-a,a, - 



3 1 






(( a 2 +a) 2 a 3 1 /2 



j 

0 , o 1|2 0 2 -1/2, s a 0 x 

-3 (a 2 +a) aj a^So^a! (a 2 +a)~ — } 

a. 



(142) 



r 21 = 



a o a i 



3/2 



2 3 

a 0 + a l 



1 «1 



a. 



• 2 *L a ' 



3^81, -o — i 



6a, 



(143) 



a. 



r 22 = 



l2*2 2 3 

* a 0 + a l 



(a 2 + a) - — 
a i. 



(144) 



r 23 “ 



1 « 



5/2 



,2.2 2 3 

ft 6^ a 0 + ai 



(-(a 2 +a) 2 ai~ao+2a 1 a 0 (a 2 +a)) 



(145) 



r 24 = 



a o a i 



3/2 5/2 



■[a 



a. 



1 «, 



al + ccf a o 3b 2 8 ^ a o 



(146) 



/, s3 3/2 N 2 1/2 o 2 -1/2, s a 0 i, 

t(« 2 + «) «i -3( a 2 +aYa 1 a 0 + 3a 0 a 1 (a 2 + a) — )] 



a, 



The cubic coefficients K and M are 



K = 0.03603 [ -0.733 1 + 6.7779 a 3 

+—!— (6.7779 cd* + + 14.5725 <d^ ( -cd^ + (4.3 cdJ - 1 ) ( 1.75 cd„ + a)) 

ft 2 6^ 1-4662 cd„ 

+ 07166 (1.884 cd* + 1.4663a to) 3 )] 

(147) 
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M =0.036031--^ + 6.45a, (o*3 
4.3 3 " 

+ — — [ -21.3675o) 3 +2.15(0 3 ( -7.6314(0* - 1 1.8787a co +2.15 2 a 2 (o^ 148) 

Z2, 4 



- 2. 15 2 (o fi ( 1 .884 oo n + 2. 1 5 1/2 a (o n ) 



1 



4.3988(0 



■(1.884(o 3 + 2.15 1/2 ao) /1 ) 3 ]] 



5. Perturbation in b 



The equations of motion are, 



vjr = r 
r = a 
y = sinty 



r = ar+a 3 r 3 +cbb 



(149) 



where, 



and the Jacobian is, 



6 = 6_„,tanh 

sat 



M 

°sat) 



5 0 = +K r r+K y y 



(150) 



A = 



0 1 0 
-cdj (l-c)a-co 2 -ca 0 

1 0 0 



(151) 



The eigenvalues are, 
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(152) 



x ia = 



a 0 +aOj 

>1 a 2 +a 



X 3 - 



a r 



a. 



The normal coordinate transformation is, 

= z x 



r = - 



\ 



a +a a. 

o 1 



ct 2 +a 



y = 



a 2 +C 

\| a o +aa i 



The normalized equations of motion are, 



ij = p'ez,- 



/ 

/ 


\ 

a 0 +aa, 


CO G + 




N 


a, +<t 
2 / 



o e + 



*2 = 

where /?’ and «' are, 



a 0 +aa l 



N a 2 



+a 



h +r n z l +r n z \h + r n z \ z l +r u z l 



z { + P'eZj +r 21 z 1 3 +r 22 z 1 2 z 2 +r 23 z 1 z|+r 24 z 2 3 



P' 



ai(a 0 + aja)(a 2 +a) 

2(a 0 + a i a ) a i + a o( a 2 +a ) 



2.15 2 (Q n +2.15a)(1.75w n +a) 
10.9950 +20.8767a 

n 



(153) 



(154) 



(155) 
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a l 


a 2 +a 


a 0 ( a o +aa i)( a 2 +a ) ' 

X 




a 0 +aa, 


v (a 0 +aa l )a 2 1 + a 2 0 (a 2 + a) y 



= 1.075 



\ 



1.75 co n + a 
o) n +2.15 a 



^ (fa) n +2.15a)(1.75fa) /i +g) \ 

(6.3725 + 10.9384a) j 



(156) 



The coefficients of third order terms are given in terms of the dummy coefficient fl , 



“i« 0 (ct 2 +a) 1 / 2 (a 0 + aa 1 ) 1/2 

aj(a 2 +a)+ai(a 0 +aa 1 ) 



(157) 



7*n = Q 



1 

6 N 



a 0 +aa x a x 



*2 +a 3b 2 b 2 sat \ 



a 0 +aa l 
a 2 +a 



r\ a 3 

r n = C 1 —TT a i 
b 2 b 2 



'sat 



7*,c =0 



a 2 a 2 



1 ^ 2 6LN 



«2 +g 

a 0 +aa 3 



r !4 = Q 






a„ + aa 



a 2+ a f 



(aa^ 



3/ 



3b 2 b 



2* 2 



a n +<za 



\ 3 

i 



SOf V 



a 2 



r 21 = " Q 



a o 1 “1 


r a 0 +aaj\ 


6a i 3b 2 b 2 sat a o 


a, +a 
V 2 / 



(158) 



(159) 



(160) 



(161) 



(162) 
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r 72 =Q 



C cru a «1 

b 2 b L 



a Q +fla 1 
a 2 + a 



2 4 



r 23 =" Q 



a a 



2> 2 6 



sat 



(163) 



(164) 



r 24 = Q 



a. 



a, 



a. 



a, +aa 



\3/2 



a 2 + a 



a 



3 a. 



3b 2 b 2 sai a 2 0 \ 



“ 2 +a 
a 0 + fla, 



(165) 



The cubic coefficients K and M are expressed in terms of the natural frequency o, 



K = 



o„ 



8(6.3725 o> 2 + a oo n + 2. 15 2 a ) 
2.15 3 (o 2 



[2.15a J o n 



(<o +2.15a) 2 

J G) 

1.75a> + a 



-1.075(0+2.15 a) 



(166) 






- [2.15o)„(o n +2. 15a) 2 +a(1.75 a)/a)(l + 2.15 2 a 2 co„)J] 



M = 



0.2688 «/(l. 75 a> n +a)(a> n + 2.15a) <o n +a2.15 co 

- [3a, 



(2.9 co + 6.6225a) 



1.75 o) + a 4.3 

tt 



b 2 b 



_ . 3 <i> +2.15a 

-I2.15 4 o)^ — 

1.75o) +a 



-(2.15a 4 <o* 



SOT 



-2.15 3 acof +2.15 3 a 3 coJ 2 



a> +a2.15 \ 3 



1.75 w n +a 



(167) 
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C. RESULTS 



Before we present the results of analytic study of Hopf bifurcation we have to 
assess the effect of the a 3 term. The a 3 term plays a role in equation (12), and from the 
nature of the equation we can see that a 3 has to have a negative value. This is because 
of the softening spring characteristic of the steady state r - 6 curve, as shown in 
Figure 10. 




5 

Figure 10 r -6 graph for different a 3 values. 



In the previous sections expressions for /S’ and K were found for equation (89). 
The value of K depends on only the nonlinear terms and since the eigenvalues of the 
Jacobian matrix cross the imaginary axis with nonzero speed the term )3 ' is nonzero. 
Equation (89) has two steady state solutions, one at R=0 which corresponds to the 
trivial equilibrium solution at zero and one at 
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(168) 




This equilibrium solution corresponds to a periodic solution (limit cycle) in the cartesian 
coordinates z,, From equation (168) we can conclude that, 

1 . If 6 ' > 0 then, 

a. If K > 0, then unstable period solutions coexist with the stable equilibrium 
for e < 0, and 

b. If K < 0 , then stable period solutions coexist with the unstable equilibrium 

for e>0, 



2. If 6 ' <0 then, 

a. If K>0, then unstable period solutions coexist with the stable equilibrium 
for e>0, and 

b. If K < 0 , then stable period solutions coexist with the unstable equilibrium 

for € < 0. 

The stable periodic solutions form the supercritical Poincare ' - Andronov - Hopf 
(PAH) bifurcation, while the unstable periodic solutions form the subcritical PAH 
bifurcation. The period of the limit cycles is computed by substituting equation (168) in 
(94) 



rp 2 tc 2n to'K-ft'm . 2 \ 

T= = — (1 — e) +0(e z ) 



o) n + oi'e +MRq “n 






( 169 ) 
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c 



A 




Figure 11 The two general types of Hopf bifurcation. 



The existence and stability of periodic solutions in our cases is examined in the 
following subsections. The graphs of u' and 0' are presented in figures (12) and (13). 
We note that for b and for natural frequency in the range | a/1 .75 1 <oo n < 1 2. 15a | the 
system does not have a complex pair of eigenvalues. Hopf bifurcation does not occur in 
that interval. 



46 




Wn 



Figure 12 /? ' versus co n . 




Figure 13 w' versus cj q . 
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1. Perturbation in 

The graph of fi ' for is shown in figure (13). /S' is always negative for K^, 

and the solutions of equation (91) is presented in figures (14), (15), (16). K is less than 
zero for low natural frequencies. For those frequencies we get a supercritical Hopf 
bifurcation and a stable periodic solution exists. However when the natural frequency 
increases the bifurcation shifts to a subcritical Hopf bifurcation. Figure (14) is obtained 
for a 3 = 0 and various o^,. Figure (15) is obtained for a 3 = -3 and various 5 sat . Figure 
(16) is obtained for o sat =0.4 and various a 3 .When |a 3 | increases the domain of the 
natural frequencies, which causes supercritical Hopf bifurcation, increases as well. 




Wn 

Figure 14 K Ki versus « n for a 3 = 0 and various S sat . 
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Wn 

Figure 16 K k ^ versus oo n for 5 sat — 0.4 and various a 3 . 
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2. Perturbation in K r 

The term $' for K r changes its sign at w n = |a/1.75| and can be observed 
from Figure (12). Therefore, for values of K<0 there exists a supercritical Hopf 
bifurcation whereas subcritical Hopf bifurcation forms when K changes its sign. From 
Figure (17) it is observed that 5 Mt does not affect the domain of the natural frequency 
for supercritical Hopf bifurcation. However, an increase in the |a 3 | term increases the 
domain of supercritical Hopf bifurcation, as can be seen from Figures (18) and (19). 




Wn 

Figure 17 K Kr versus co n for a 3 =0 and various S sat . 



50 





Wn 

Figure 19 versus a) n for 6 sat =0.4 and various a 3 . 
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3. Perturbation in Ky 

For this case /3 ' is always positive and for e > 0 there exists an unstable 
equilibrium, therefore for K<0 a supercritical Hopf bifurcation exists. Similarly for 
K > 0 a subcritical Hopf bifurcation exists. An increase in 5 sat increases the domain of the 
supercritical Hopf bifurcation over the natural frequency, see Figure (20). Also an 
increase in | a 3 1 increases the domain of the supercritical Hopf bifurcation, see Figures 
(21), (22). 




Figure 20 K Ky versus u) n for a 3 = 0 and various 5 sat 
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Figure 21 K Ky versus co n for a 3 = -3 and various 8 , 



sat 
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4. Perturbation in a 



( 3 ' is always positive for this case. For e > 0 and K < 0 there exists a stable 
period solution coexisting with the unstable equilibrium. When a 3 = 0, a change in 5 sat 
does not affect the supercritical Hopf bifurcation region of the natural frequency, Figure 
(23). If a 3 is non zero, an increase in 8 ^ increases the domain of the natural frequency 
which causes the supercritical Hopf bifurcation. An increase in |a 3 | also increases the 
domain of the natural frequencies for supercritical Hopf bifurcation, see Figures (24), 
(25). 
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Figure 24 K a versus a> n for a 3 = -3 and various 5 sat 




Wn 

Figure 25 K a versus w n for o sat — 0.4 and various a 3 . 
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5. Perturbation in b 



Previously it was shown that there is no Hopf bifurcation in the interval, 
1 1.75/a | <oj n < 1 2. 15a | K can be examined in two regions, 0<co n < 1 1.75/a | and 
w n > 1 21 . 15a | . In the first region ' is always negative, e > 0, and K is always negative. 
In this region only supercritical Hopf bifurcation occurs. Changes in 6 sat and |a 3 | have 
no effect in this region, see Figures (26), (27), (28) . For the second region (3 ' <0, e <0, 
K<0 initially and supercritical Hopf bifurcations exist. For higher natural frequencies, 
subcritical Hopf bifurcation is observed. Figures (29) and (30) show that changes in 5 sat 
or a 3 have no significant effect. 




wn 

Figure 26 K b versus w n for 5 sa , =0.4 and various a 3 for the first region. 



56 




Figure 27 K b versus oj n for a 3 = 0 and various 5^ for the first region. 




Figure 28 K b versus co D for a 3 = -3 and various <5 sal for the first region. 
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wn 

Figure 30 K b versus co D for 5 sat = 0.4 and various a 3 for the second region. 
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IV. SIMULATIONS 



To illustrate the dynamics of the Hopf bifurcation we simulated the case in which 
there is a change in the gain Ky. We compared the simulations with the analytic studies 
presented in previous chapters. In our case, loss of stability occurs at c=3.7625 from 
stable to unstable while c increases its value, as can be observed from Figure (6). In the 
following simulations we use the typical saturation angle for the rudder of 0.4 radians. 
The a 3 term is assumed to be zero. 

K is less than zero for the natural frequency 0.4 in Figure (20), therefore a 
supercritical Hopf bifurcation exists in this region. The behavior of this supercritical 
Hopf bifurcation is simulated in Figures (31), and (32). In Figure (31) c Ky =3.5 is 
simulated for two different initial lateral deviations ( y 0 = 0.05 and y 0 = 0.5). 
Regardless of the initial conditions both curves converge to zero for this value of c. 
There is a unique stable steady state solution for this value of c. Figure (32) is simulated 
for K < 0 and a> n = 0.4 and the same initial condition of y 0 =0.5. When c= 3.5 the 
simulation converges to zero as explained in Figure (31). When c is increased above the 
critical value (here c=4.0) we observe that the solution converges to a limit cycle. This 
stable limit cycle coexists with the unstable equilibrium solution. 

Figures (33) and (34) are presented to show the effects of a subcritical Hopf 
bifurcation. For this purpose w n =1.2 is simulated for Figures (33), and (34). u n = 1.2 
is in the region where K>0 in Figure (20). The initial condition for this case is 
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y 0 = 0.05. When c is set to 3.5 the solution converges to zero, but for the case where 
c is set to 4.0 the simulations converge to a stable limit cycle. This is the case in Figure 
(11-b). The simulation misses the unstable limit cycles and converges to the outer limit 
cycle shown in Figure (11-b). In Figure (34) the effects of the initial conditions are 
further illustrated. In Figure (34) for y 0 =0.05, the simulations converge to zero, but for 
y 0 =0.5 the simulations converge to a limit cycle. For this simulation the initial condition 
was large enough to shift the results to a limit cycle rather than the zero steady state 
solution. 

From the simulations we observe that in subcritical bifurcations the magnitude of 
the stable limit cycle is greater than the magnitude of the limit cycle in supercritical 
bifurcations. Similar results can be obtained for the other cases of variation in parameters 
a, b, and gains and K y . 
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SUPERCRITICAL HOPF BIFURCATION 




Figure 31 Supercritical Hopf bifurcation for oj n =0.4 , c=3.5 for two simulations 

with y 0 =0.05 and y 0 =0.5. 



SUPERCRITICAL HOPF BIFURCATION 




Figure 32 Supercritical Hopf bifurcation for cj n =0.4 , y 0 = 0.5 for two simulations 
c=3.5 (stable equilibrium ) and c=4 (unstable equilibrium) . 
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SUBCRITICAL HOPF BIFURCATION 




Figure 33 Subcritical Hopf bifurcation for oj n = 1 . 2 , y 0 = 0.5 for two simulations 
c=3.5 (stable equilibrium ) and c=4 (unstable equilibrium) . 



SUBCRITICAL HOPF BIFURCATION 




Figure 34 Subcritical Hopf bifurcation for co n = 1.2 (K>0) where c= 3. 5. The simulations 
for y o =0.05 converged to the steady state value and for y 0 =0.5 to the limit cycle. 



62 



V. SUMMARY AND CONCLUSIONS 



An analytic investigation of the nonlinear dynamic response characteristics of a 
steering control law of marine vehicles has been presented. Bifurcation theory techniques 
were utilized in order to assess the behavior of the system upon initial loss of stability. 
The main bifurcation parameters were the natural frequency, control saturation level, 
system parameters and gains. The main conclusions of this research can be summarized 
as follows: 

1. There exists a critical point for each gain and system parameters for stability of 
straight line motion. The loss of stability occurs in the form of Hopf bifurcation for each 
case. As the parameter crosses its critical value, a family of periodic orbits ("self 
oscillations") develops. 

2. The desired characteristic equation from the ITAE criterion has a faster response 
for higher natural frequency. However, subcritical Hopf bifurcation develops when the 
natural frequency is sufficiently high. If c is increased past c crit the amplitude of lateral 
deviation undergoes a jump and takes "large" values ( "hard loss of stability or hard 
generation of limit cycles"). For lower natural frequencies, stability loss occurs as 
supercritical where the limit cycle amplitude grows continuously (" soft loss of 
stability"). 

3. The stationary point where change of stability occurs is related only to the natural 
frequency of the system and the hydrodynamic coefficient a. 
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4. Higher saturation limit of the rudder angle increases the range for supercritical 
Hopf bifurcation. Theoretically, if a saturation limit does not exist all bifurcations 
develop as supercritical. 

5. The term a 3 has a similar effect as the saturation limit of rudder angle. An increase 
in the |a 3 | term increases the natural frequency where supercritical Hopf bifurcation 
exists, but an increase in the | a 3 1 term is not a desired condition since vehicle response 
to rudder angle slows down. 

6. Transitions from supercritical to subcritical are very important for the design of the 
control law. A supercritical Hopf bifurcation is preferred over a subcritical Hopf 
bifurcation since the subcritical Hopf bifurcation may develop a rapid dynamic jump to 
a new limit cycle. 

Some recommendations for further research are as follows: 

1. Comparative studies must be performed including the observer dynamics. The 
effect of variations in gains and system parameters must be studied to predict the 
accuracy of the observer. 

2. Continuation techniques for periodic solutions must be performed for further 
stability analysis. The cases for simultaneous variations of more than one terms must be 
studied to understand the dynamic response of the system including any nonlinear terms. 
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